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Abstract 

We consider two models of interacting Bose gases: a gas of spin one particles in the ground 
state of a cubic box and a one-dimension Bose gas with contact interactions. We show how 
to calculate exact eigenstates of the corresponding N— body Hamiltonians. Both models 
share the property of not leading to the formation of a Bose-Einstein condensate, even at 
zero temperature, in the strict sense of the existence of a single one-particle state with a 
macroscopic population. We show that a lot of physical insight can be gained on these two 
model systems by using the usual Hartree-Fock mean field approach: in this approximation, 
that we test against the exact result, everything happens as if a single realization of the 
system was a Bose-Einstein condensate in a state <f> breaking the rotational or translational 
symmetry, and varying in a random way for any new experimental realization. 

Since the observation in 1995 of the first Bose-Einstein condensates in atomic gases ^ |2J E] 
a renewed interest is taking place for these macroscopic quantum states of matter . 

In particular some experiments are the realization of what was considered before 1995 as 
gedanken experiments. From this point of view the interference experiment performed at MIT 
by the group of Wolfgang Ketterle is particularly illustrative: two condensates that 'had never 
seen each other' are made to spatially overlap, leading to the formation of interference fringes jSJ. 
In this last sentence we have deliberately paraphrased P.W. Anderson, who raised the question 
whether two superfluids that 'have never seen each other' have or not a well-defined relative 
phase 

The MIT-type phase experiment and more generally the concept of phase of a matter wave 
field has been the subject of several theoretical analyses (ZHHUHlEl- One of the theoretical issues 
was to reconcile two different points of view. In the first point of view, the two condensates 
in the MIT experiment are described by a Fock state, that is by a state with a well-defined 
number of particles in each condensate mode; this Fock state leads to interference fringes in a 
measurement of the positions of all the atoms, as shown numerically in [Hj and analytically in 
|1U| . but the result is not obvious. In the second point of view a symmetry breaking descrip- 
tion is adopted, attributing a well-defined relative phase cj) to the two condensates in a given 
experimental realization, with <p varying in an unpredictable manner for any new realization 
of the experiment, with a uniform probability distribution in [0, 2tt[ . The symmetry breaking 
point of view immediately predicts the formation of interference fringes, in a very intuitive way, 
but requires some justification, as a symmetry of the system is broken! We refer the reader to 
|10[ II 1 j for a comparison of the two points of view, the main result being that their predictions 
for the interference fringes coincide in the limit of large number of particles, and that symmetry 
breaking is a convenient but by no means necessary procedure. 

* Laboratoire Kastler-Brossel is a unite de recherche de l'Ecole normale superieure et de l'Universite Pierre et 
Marie Curie, associee au CNRS. 
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The goal of the present text is to illustrate other uses of symmetry breaking, allowing one to 
extend the applicability of mean-field approximations to situations where there is no condensate 
in the strict sense, that is in the sense of Penrose and Onsager |12j . but where a so-called 
fragmented condensate ^3] is present (section The symmetry-breaking point of view that 
we put forward here is of a physically different nature than the U(l) phase symmetry breaking 
description of the above mentioned MIT interference experiment. We consider in section [21 the 
case of so-called 'spinor condensates', more specifically condensates of particles of spin one, 
where a SO (3) rotational symmetry breaking description is applied. We treat in section 01 the 
case of one-dimension attractive Bose gases, where a spatial translational symmetry breaking is 
applied to obtain solitonic condensates. 

In both cases of spinor and solitonic condensates the procedure followed in this paper is 
the same one: the ground state of the system is calculated exactly for an appropriate model of 
the interaction potential and it is of course symmetric, while its mean-field approximation by 
Hartree-Fock states breaks the symmetry. In both cases we will consider Gedanken experiments 
whose single outcomes can be predicted easily from the Hartree-Fock state and with more effort 
from the exact ground state, the two predictions coinciding for a large number of particles. This 
will illustrate the ability of the mean-field approximation to allow physical predictions in an 
easy and transparent way, correct in the limit of a macroscopic system. 

1 Physical context and motivation 

1.1 Usual definition of a condensate 

We consider a gas of N interacting indistinguishable bosonic particles. The total number of 
particles is conserved by the Hamiltonian and the N— body state of the gas is defined by the 
N— body density operator > supposed to be known. When does one say that the gas is 

in a Bose-Einstein condensed state ? We answer this question in the way formalized by Penrose 
and Onsager [12] . 

One first defines the one-body density operator p\ as the trace of oi,...^ over the state of 
the particles 2, . . . , N : 

pi = iVTr 2j ...,jv [<ti,...,n] ■ (1) 

Note the factor N in this definition, so that {r\p\\f) is exactly the mean spatial density in r. 
The one-body density operator can also be defined by its matrix element in any single-particle 
orthonormal basis ip a using second quantization: 

{^ a \PMp) = (al&a) (2) 

where a Qj/ g annihilates a particle in the state ip a p and (...) stands for the expectation value 
in the density operator (Ti,...,^ • 

Second one diagonalizes the hermitian operator pi : 

Pi = ^n k \<t> k ){<t>k\ (3) 
k 

where the <pk 's form the orthonormal set of eigenvectors of p\ with eigenvalues rik ■ All 
eigenvalues are greater than or equal to zero and their sum is equal to N . One then says 
that a condensate is present if there exists a level index k , taken by convention to be k = , 
such that the two following conditions are satisfied: 



n o ^> n k f° r & H k 7^ 
no ~ N 



(4) 
(5) 
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Physically this means that the maximally occupied single particle level <t>o has a mean occupa- 
tion number uq much larger than any orthogonal single particle state (according to condition 
(H|)) and that no is macroscopic (according to condition (JSJ)), say at least on the order of 10% 
of the total number of particles in the system. 

In this case 4>${r) = {r\(fro) is the condensate wavefunction, normalized to unity, and no 
is the mean number of particles in the condensate, usually noted as No . 

1.2 Situations considered in this paper 

We consider systems such that even at zero temperature, that is even in the ground state, one 
of the two conditions (|4I5|) is violated, so that there is no Bose-Einstein condensate in the strict 
sense of [T2j . 

More specifically we will exemplify a class I of such systems corresponding to a p\ having 
several eigenstates with a macroscopic population, and a class II where none of the eigenstates 
of p\ has a macroscopic population. 

Our example for class I is a model for a gas of interacting spin one particles, all in the 
ground state of a box with periodic boundary conditions but with free spin variables. One then 
finds that because of the SO (3) rotational symmetry of the Hamiltonian the one-body density 
operator pi contains three eigenstates with the same macroscopic population (equal to N/3 ): 
condition @ is violated. 

Our example for class II is a one-dimension Bose gas with attractive contact interactions in 
free space (that is in the absence of confining external potential). One then finds that because 
of spatial translational symmetry the eigenstates of p\ are plane waves. If one normalizes these 
plane waves in a fictitious quantization box of size L one finds that all plane waves have a 
population tending to zero in the limit L — ► +oo : condition (JSJ) is violated. 

1.3 Methodology of the symmetry breaking point of view 

In the usual case of a Bose-Einstein condensed system in the Penrose and Onsager sense, p\ 
has a single macroscopic eigenvalue No with a corresponding normalized eigenvector 1 0o) • 
When the gas is dilute, that is when the mean interparticle separation is much larger than the 
scattering length a of the interaction potential, the fraction of particles out of the condensate 
1 — Nq/N is very small at zero temperature and the interaction potential can be replaced by 
an effective low energy contact interaction |1J . One can then determine an approximation 
to <f) by using the following Hartree-Fock "pure condensate" variational ansatz for the ground 
state of the gas: 



that is a Fock state with all the N particles in the same mode (ft . One determines <p by 
minimizing the mean energy of . This leads in particular to the famous Gross-Pitaevskii 

equation for 4> > verv successful in predicting the properties of the Bose-Einstein condensates 
of alkali gases [4-. Note that after minimization there is still the freedom of a phase change of 
4> by a constant phase factor e ia ; but this amounts to multiplying the N— body state vector 
|^ r ) hf by a global phase e lNa of no physical consequence. 

We shall use exactly the same Hartree-Fock variational procedure in our examples of a spin 
one Bose gas and of a one-dimension Bose gas with attractive interactions. Something peculiar 
will happen then: one finds an infinite number of minimal energy solutions, rather than a single 
one (up to a global phase factor) in the usual case. These minimal energy solutions have all 
exactly the same energy, and are parameterized by a set of continuous parameter(s) Vt . We 
note (fin the minimal energy solution corresponding to the parameter(s) . In the case of spin 
one particles is actually a particular direction in space in which the spin of </>n points. In 
the one-dimension Bose gas £1 gives the position of the center of mass associated to </>q . 




(6) 
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This structure of the results is a direct consequence of the fact that minimal energy solutions 
of the Hartree-Fock ansatz are not symmetric with respect to the symmetry group of the Hamil- 
tonian, as they privilege some arbitrary spin direction (for the spin one gas) or are a localized 
soliton around some arbitrary point of space (for the one-dimension Bose gas). Each cpQ is 
said to break the symmetry, the rotational SO (3) symmetry for the spins or the translational 
symmetry for the ID Bose gas. From any given 4>q one then obtains a continuous family of 
(j)Q, with the same mean energy by applying the continuous symmetry group of the Hamiltonian 
to (fttt j that is arbitrary spin rotations or arbitrary spatial translations. 

How can we then construct a reasonable approximation to the exact ground state |^o) °f 
the gas from this continuous family of Hartree-Fock states? The simplest choice is the so-called 
symmetry-breaking prescription, approximating the pure state |^o) by the following N— body 
density operator: 

»t,N = J ^V:0n)(iV:M (7) 

that is by a statistical mixture of all the minimal energy Hartree-Fock ansatz. The normal- 
ization factor V is independent of O for a convenient choice of the parameterization of the 
symmetry group [Tl]: e.g. for SO (3) one integrates over all possible solid angles . Why is the 
statistical mixture (JJJ) called a symmetry breaking description, while it is obviously invariant 
by the symmetry group ? Well, this prescription allows one to imagine that each particular 
experimental realization of the gas corresponds to a condensate with all particles in the same 
state 4>q, , where varies in a random, unpredictable way for any new realization of the ex- 
periment: each particular experimental realization therefore 'breaks' the symmetry. This point 
of view should of course not be taken too literally: it is only a convenient reinterpretation of the 
density operator 0. The same density operator can also be obtained as a statistical mixture 
of other N— body states than the (f>n (e.g. a statistical mixture of its eigenstates) which may 
lead to a very different physical picture. What matters actually is that the quantum mechanical 
prediction of the results of any measurement performed on the gas will be expressed as the 
expectation value of some operator in the density operator d"i jv |1U1 111 j : 

(6) = Tr [dai,...,jv] (8) 

so that physical predictions depend only on &i n as a whole, not on a particular decomposition 
of o"i jv as a statistical mixture. For example the probability density of measuring the particle 
1 in point of real space r± , . . . , the particle N in point r/v corresponds to the expectation 
value of the N— body observable 6 = &(r N ) . . . *t(ri)#(ri) . . . &(r N )/N\ , where * is the 
atomic field operator |Hj. 

There exists a more elaborate way to approximate the N— body ground state of the gas by 
a pure state constructed from the Hartree-Fock states \N : <po) . The idea is to form a linear 
combination of these Hartree-Fock states having the same symmetry as the exact ground state. 
This is achieved with the use of a projection theorem associated to the representations in the 
Hilbert space of the symmetry group of the Hamiltonian For a continuous symmetry group 
this amounts to performing integrals of the form 

vhfs 



|* r s = J dQf(Q)\N : <h) (9) 

for an appropriate choice of the weight factor /(O) . We restrict for simplicity to the case 
of a non-degenerate ground state, associated to a unidimensional irreducible representation of 
the symmetry group; the weight factor /(O) then has a modulus independent of £1 . While 
mathematically more satisfactory than this form is not particularly illuminating from a 
physical point of view. It is also more difficult to manipulate. Let us indeed try to calculate the 
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expectation value of an arbitrary k -body operator O in the state ©. In first quantized form 
O is a sum over all possible k -tuplets of particles of an operator of k particles: 

6= Yl o(h,...,i k ). (io) 

l<ii<...<i k <N 

( N 

Using the bosonic symmetry of the Hartree-Fock states one finds that all the I , 
N\/[k\(N — k)l] possible k -uplets give the same contribution to the expectation value so that 

-<*o|o|*or =( N k )f <mm) I ^ r( n <fc^W--- ; fc)lfc^> m ,\ M f. 

V k ) J J ((<ki>\<t>n)) 

(ii) 

In general the double integral over O and O' in Eq. ([ll|) is difficult to calculate. In the case 
k <C N , one can fortunately use a large iV expansion as follows. The scalar product of 4>q 
and 4>Q' raised to the very large power N is a much narrower function of the distance £1 — 0,' 
than the remaining part of the integrand so that it can be approximated by a 5 distribution 
QUI: 

((^1^)^-^(0-0') (12) 
where J\f is a constant factor. The expectation value (|TT|) then reduces to 

Ms (* |O|* ) Ms ~ ( N k ) N j dn\f(n)\ 2 (k:cfi n \0(l,...,k)\k:^}. (13) 

As | /| has a constant modulus, this expression is equivalent to the symmetry-breaking predic- 
tion 0! 

One could then go to next order in the 1/N expansion, to identify the differences in the 
predictions of ((7J) and ©. In general and as we shall see for the one-dimension Bose gas, this is 
not worth the effort as both hfs ( 1 I / o|0| v I / o) llfs and the symmetry-breaking prediction Tr [Ocx^ N ] 
differ from the exact result by a term of the same order of magnitude, of order 1 /N times the 
exact result. An amusing exception is our model of spin one particles where © actually coincides 
with the exact ground spin state! 



2 The ground state of spinor condensates 

The alkali atoms used in the Bose-Einstein condensation experiments have an hyperfine structure 
in the ground state, each hyperfine level having several Zeeman sublevels. Consider for example 
23 Na atoms used at MIT in the group of Wolfgang Ketterle. The ground state has an hyperfine 
splitting between the lower multiplicity of angular momentum F = 1 and the higher multiplicity 
of angular momentum F = 2 . All the three Zeeman sublevels mp = 0, ±1 of the lower 
multiplicity F = 1 cannot be trapped in a magnetic trap (if rap = — 1 is trapped then 
mp = +1 which experiences an opposite Zeeman shift is antitrapped). But they can all be 
trapped in an optical dipole trap, produced with a far off-resonance laser beam, as the Zeeman 
sublevels experience then all the same lightshift. This optical trapping was performed at MIT 
|16j . opening the way to a series of interesting experiments with condensates of particles of spin 
one P2|. 

We introduce first a model interaction potential for spin one particles. We then concentrate 
on a specific aspect, the ground state of the spinor part of the problem, assuming that the atoms 
are all in the ground state of a cubic box with periodic boundary conditions. 
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2.1 A model interaction potential 



In the theory of Bose-Einstein condensates the interaction potential between the particles is usu- 
ally replaced by a model interaction potential having the same collisional properties at low en- 
ergy as the 'exact' atomic interaction potential. This model interaction potential is conveniently 
taken to be Fermi's pseudo-potential initially introduced in nuclear physics [181 1191 QUI lllj . It 
has the following action on a two-body wavefunction ipi 2 : 



(ri,r 2 |F(l,2)|^ li2 ) = g8{f{-r 2 ) 







(V12 ^1,2(^1,^*2) 



(14) 



This is essentially a contact interaction potential, conveniently regularized, and with a coupling 
constant g proportional to the s -wave scattering length a a of the 'exact' interaction potential: 

.2 
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-a. 



m 



(15) 



We have to generalize the model scalar pseudo-potential of Eq. (|14|) to the case of particles 
having a spin different from zero. As we want to keep the simplicity of a contact interaction 
potential we choose the simple form 



V(l,2) = Vspin(l s 2)<5(r!-r5) 



d 



Or 



12 



(ri2 



(16) 



that is the product of an operator acting only on the spin of the particles 1 and 2, and of the 
usual regularized contact interaction acting only on the relative motion of the two particles. 
The interaction potential V(l,2) has to be invariant by a simultaneous rotation of the spin 
variables and of the position variables of the two particles. As the contact interaction is already 
rotationally invariant, the spin part of the interaction V sp i n (l,2) has to be invariant by any 
simultaneous rotation of the two spins. 

This condition of rotational invariance of V sp m(l,2) is easy to express in the coupled basis 
obtained by the composition of the two spins of particle 1 and particle 2: within each subspace 
of well defined total angular momentum V sp i n (l,2) has to be a scalar. Let us restrict to the 
case studied at MIT, with spin one particles. By composition of F = 1 and F = 1 we obtain 
a total angular momentum F tot = 2 , 1 or 0, so that one can write 



V spin (l, 2) = g 2 P Ft , =2 (l, 2) + 5l P Ftot=1 (l, 2) + g P Ftot =o(l, 2) 



(17) 



where the g 's are coupling constants and the P(l,2) 's are projectors on the subspace of 
particles 1 and 2 with a well defined total angular momentum F to t . At this stage we can play 
a little trick, using the fact that the states of Ftot = 1 are antisymmetric by the exchange 
of particles 1 and 2 (whereas the other subspaces are symmetric). The regularized contact 
interaction scatters only in the s -wave, where the external wavefunction of atoms 1 and 2 is 
even by the exchange of the positions r\ and r 2 ; as our atoms are bosons, the spin part has 
also to be symmetric by exchange of the spins of atoms 1 and 2 so that the 'fermionic' part of 
V sp in(l, 2) , that is in the subspace Ftot = 1 , has no effect. We can therefore change g± at will 
without affecting the interactions between bosons. The most convenient choice is to set g\ = g 2 
so that we obtain 

V spin (l, 2) = 52 Id(l, 2) + ( 50 - g 2 )P Ftot =o(l, 2) (18) 

where Id is the identity. The subspace Ftot = is actually of dimension one, and it is 
spanned by the vanishing total angular momentum state \ipo(l,2)) . Using the standard basis 
|m = —1,0,-1-1) of single particle angular momentum with z as quantization axis, one can 
write 

hfo(l,2)> = — [1 + 1,-1) + I -1,+1)-|0,0)]. (19) 
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A more symmetric writing is obtained in the single particle angular momentum basis \x,y,z) 
used in chemistry, defined by 

1 + 1) = -±={\x) + i\y)) (20) 

1-1) = +^(\x)-i\y}) (21) 
|0) = \z). (22) 

The vector \a) in this basis (a = x,y,z) is an eigenvector of angular momentum along axis 
a with the eigenvalue zero. One then obtains 

|Vo(l, 2)) = [\x, x) + \y, y) + \z, z)} . (23) 

To summarize, the part of the Hamiltonian describing the interactions between the particles 
can be written, if one forgets for simplicity the regularizing operator in the pseudo-potential: 



Hint = y J d 3 f ^l^l^^a 
~ a,(3=x,y,z 

+ 90_92j d 3 f J2 fafaMfi (24) 

a,/3=x,y,z 

where ip a (r) is the atomic field operator for the spin state \a) . This model Hamiltonian has 
also been proposed by |2T1 [2"2"1 [2"3"j . 

We now restrict to the spinor part of the N— body problem, assuming that all the N 
particles of the gas are Bose-Einstein condensed in the ground state of a cubic box of size L 
with periodic boundary conditions. This is a reasonable assumption for a dilute Bose gas with 
effective repulsive interactions ( g<i , go > ) . It amounts to approximating the atomic field 
operators by 

where the operator a a annihilates a particle in the state \k = 0)\a) (a = x,y,z) that is in 
the plane wave of vanishing wavevector in the box and with internal spin in state \a) . This 
leads to our final model Hamiltonian: 

#s P in = ^3 2 a) a a\a p a a + ^ (g - g 2 )A^ A (26) 

a,/3=x,y,z 

where we have introduced 

A = a 2 x + a 2 y + a 2 z . (27) 

Up to a numerical factor A annihilates a pair of particles in the two-particle spin state |^o(l>2)) 
of vanishing total angular momentum, as shown by Eq. (|2«-Sj) . Note that we suppose that there is 
no magnetic field applied to the sample; the effect of a magnetic field is considered in |24| . 



2.2 Ground state in the Hartree-Fock approximation 

We now minimize the energy of the iV spins within the Hartree-Fock trial state vectors | N : <f>) 
with the constraint that \<f>) is normalized to unity but without any constraint on the total 
angular momentum of the spins. The external part of the condensate wavefunction is simply 
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the plane wave with momentum k = whereas the spinor part of the wavefunction remains to 
be determined: 

( f \^ = jjj2 E c » with Ei c «i 2 = L ( 28 ) 

a=x,y,z a 

From the model interaction Hamiltonian (|24|) the mean energy per particle in a Hartree-Fock 
state is easy to find, using the identity making Hartree-Fock states so convenient: 

# a (?)\N : 4>) = N 1/2 (f, a\(j))\N - 1 : <j>) (29) 

' c a \N -l:<j)), (30) 



L 3 



One then finds the mean energy per particle: 



E N-l N-l, „, 2 , s 

92 + -^rr(go -92)\A\ 2 (31) 



N 2L 3 ^ 6L 3 
where we have introduced the complex quantity 

A= 4 = ^ (32) 

a=x,y,z 

where c is the vector of components (c x ,c y ,c z ) . We have to minimize the mean energy over 
the state of the spinor. 

• Case g 2 > go 

This is the case of sodium . As the coefficient go — 92 is negative in Eq. (|3*T|) we have to 
maximize the modulus of the complex quantity A . As the modulus of a sum is less than the 
sum of the moduli we immediately get the upper bound 

\A\< £ \c a \ 2 = l (33) 

a=x,y,z 

leading to the minimal energy per particle 

E N - 1 JV-1 , 

N = ^ir 92 + -Qir {90 - 92) - (34) 

The upper bound for \A\ is reached only if all complex numbers c 2 , have the same phase 
modulo 2-7T . This means that one can write 

c a = e i6 n a (35) 

where 6 is a constant phase and n = (n x ,n y ,n z ) is any unit vector with real components. 
Physically this corresponds to a spinor wavefunction being the zero angular momentum state 
for a quantization axis pointing in the direction n . The direction n is well defined in the 
Hartree-Fock ansatz, but it is arbitrary as no spin direction is privileged by the Hamiltonian. 
We are facing symmetry breaking, here a rotational SO (3) symmetry breaking. The symmetry 
breaking prescription ® then leads to the density operator 

<..,tv = J ^\N:0 n )(N: n \, (36) 

where the integral is taken over the unit sphere, that is over all solid angles, and where we have 
introduced the single particle state 

\0h) = n x \x) + n y \y) + n z \z). (37) 
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• Case g 2 < go 

In this case we have to minimize \A\ to get the minimum of energy. The minimal value of \A\ 
is simply zero, corresponding to spin configurations such that 



c 2 - ^ - 2 



E c '=° (38) 

a=x,y,z 

with an energy per particle 

E N - 1 

N ~ Mr*- <39) 
To get more physical understanding we split the vector c as 

c = R + if (40) 

where the vectors R and / have purely real components. Expressing the fact that the real 
part and imaginary part of c 2 vanish, and using the normalization condition c • c* = 1 in 
Eq. (j28|) we finally obtain 

RI = (41) 
R 2 = I 2 = \- (42) 

This means that the complex vector c is circularly polarized with respect to the axis n or- 
thogonal to I and R . Physically this corresponds to a spinor condensate wavefunction having 
an angular momentum ±h along the axis n . The direction of axis n is well defined in the 
Hartree-Fock ansatz but it is arbitrary. The symmetry breaking prescription ([7J) then leads to 
the density operator 

&t,N= I ^\N:+ n ){N:+ n \ (43) 



47T 

where \+a) is the spin state of angular momentum +fi along axis n 



2.3 Exact ground state vs Hartree-Fock for the spinor model 

The Hamiltonian (|26|) can be diagonalized exactly |25| . This is not surprising as (i) it is ro- 
tationally invariant and (ii) the bosonic N— particle states with a well defined total angular 
momentum Sjy can be calculated: one finds that Sjy = N, N — 2, . . . , leading to degenerate 
multiplicities of H spin of degeneracy 2SV + 1 • As a consequence the one-body density operator 
pi defined in Q calculated for the ground state of the gas is spin rotationally invariant; as it 
acts on an irreducible representation of SO (3) , of spin one, it is scalar according to Schur's 
lemma |15j : 

fa = _|fc = d)(fc = 0| 0ld spin . (44) 

We are therefore in the class I of section 11.21 

To derive the exact spectrum of the spin Hamiltonian one may use in practice the following 
tricks: The double sum proportional to g 2 in Eq. (|26|) can be expressed in terms of the operator 
total number of particles iV only, 

N = Y,a j a a a . (45) 

a 

So diagonalizing -ff S pin amounts to diagonalizing A^A ! Second the total momentum operator 

S of the iV spins, defined as the sum of all the spin operators of the individual atoms in units 
of h , can be checked to satisfy the identity 



S- S + A^A = N(N + 1) 



(46) 
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so that the Hamiltonian for N particles becomes a function of S |25| : 



#spin = - 1) + ^3 (g ~ 92) 



N(N + 1) - S ■ S 



(47) 



We recall that S • S = Sn(Sn + 1) within the subspace of total spin Sn ■ 

When g2 < go the ground state of H sp i n corresponds to the multiplicity Sn = N , con- 
taining e.g. the state with all the spins in the state |+) . In this case the N— particle states 
obtained with the Hartree-Fock approximation are exact eigenstates of H sp [ n j and the symmetry 
breaking prescription (|43|) correctly gives the projector on the maximal spin multiplicity. 

When g2 > go the ground state of H sp i n corresponds to the multiplicity of minimal total 
angular momentum, Sn = 1 for N odd or Sn = for N even. In this case the Hartree- 
Fock symmetry breaking prescription (|36|) is only an approximation of the exact ground state 
of //spin • The error on the energy per particle tends to zero in the thermodynamical limit; for 
N even one finds indeed 

SE 1 . 

-^ = -^(90-92). (48) 

But what happens if one uses the pure state, symmetry restored approximation by co- 
herently summing up the Hartree-Fock ansatz over the direction ft defined in Eg. 1)35(1 ? Assume 
that iV is even; one has then to reconstruct from the Hartree-Fock ansatz a rotationally invari- 
ant state. This amounts to considering the following normalized state for the N spins: 



|*o> Ms 



VN + lj^\N:0 H ) (49) 



The state vector |^o) hfs > being non zero and having a vanishing total angular momentum, is 
equal to the exact ground state of H sp i n : 

l*o> Ms = l*o>. (50) 

The expression (|49|l can then be used as a starting point to obtain various forms of |^q) • 
If one expresses the Hartree-Fock state as the N— th power of the creation operator ^o n « 
acting on the vacuum |vac) , and if one expands this power with the usual binomial formula, 
the integral over n can be calculated explicitly term by term and one obtains: 

|tf > = AM J | vac) (51) 

where N is a normalization factor and the operator A is defined in Eq.([27|). Formula (|51|) 
indicates that |^o) is simply a 'condensate' of pairs in the state |i/>o(l,2)) . It can be used to 
expand |^o) over Fock states with a well defined number of particles in the modes m = 0, m = 
±1 , reproducing Eq.(13) of 

To be complete we mention another way of constructing the exact eigenvectors and energy 
spectrum of H spm . The idea is to diagonalize A* A using the fact that A obeys a commutation 
relation that is reminiscent of that of an annihilation operator: 

[A, i f ] = AN + 6. (52) 

In this way acts as a raising operator: acting on an eigenstate of A* A with eigenvalue A 
and iV particles, it gives an eigenstate of A^A with eigenvalue A + AN + 6 and with N + 2 
particles. One can also check from the identity (|46|) that the action of A^ does not change the 
total spin: 

[i t ,5-S ? ] = 0. (53) 
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By repeated actions of starting from the vacuum one arrives at Ea. (|51|) . creating the eigen- 
states with N even and vanishing total spin S = . By repeated actions of A^ starting from 
the eigenstates with N = 2 and total spin S = 2 (e.g. the state | + +) ) one obtains all the 
states with N even and total spin 5 = 2. More generally the eigenstate of H sp { n with total 
spin S , a spin component m = S along z and N particles is: 

/ (N-S)/2 

\\N,S,m = S) oc (A^y \S:+l) (54) 

where IS" : +1) represents S particles in the state | + 1) . From Eq. (|54jl one can generate 
the states with spin components m = S — 1, . . . , — S by repeated actions of the spin- lowering 
operator S- = S x — iS y in the usual way. We note that formula (|54"|) was derived independently 
in jH]. 

2.4 Advantage of a symmetry breaking description 

We now illustrate the physically transparent character of the symmetry breaking prescription by 
analyzing the following gedanken experiment. Imagine that we have prepared a gas of sodium 
atoms ( gi > go ) in the collective ground spin state with an even number of particles N , and 
that we let the atoms leak one by one out of the trap, in a way that does not perturb their spin. 
We then measure the spin component along z of the outgoing atoms. Suppose that we have 
performed this measurement on k atoms, with k <C N . We then raise the simple question: 
what is the probability pk that all the k detections give a vanishing angular momentum along 
z? 

Let us start with a naive reasoning based on the one-body density matrix of the condensate. 
The mean occupation numbers of the single particle spin states \m = —1) , \m = 0) and 
\m = +1) in the initial condensate are all equal to iV/3 , see Eq. (|44|) . The probability of 
detecting the first leaking atom in \m = 0) is therefore 1/3 . Naively we assume that since 
k <^ N the detections have a very weak effect on the state of the condensate and the probability 
of detecting the n— th atom ( n < k ) in the m = channel is nearly independent of the n — 1 
previous detection results. The probability for k detections in the m = channel should then 
be 

PT VC = ^- (55) 

Actually this naive reasoning is wrong (and by far) as soon as k > 2 . The first detection of 
an atom in the m = channel projects the spin state of the remaining atoms in 

|*i>=Mao|*o> (56) 

where clq annihilates an atom in spin state m = , l^o) is the exact collective spin ground 
state and A/i is a normalization factor. The probability of detecting the second atom in m = 
(knowing that the first atom was detected in m = ) is then given by 

pi (^llE^L-iSla^i*!)' 

The denominator is simply equal to iV — 1 as l^i) is a state with iV — 1 particles. Using the 
integral form (J49|) and the simple effect of an annihilation operator on a Hartree-Fock state, e.g. 

a 2 \N : Off) = [N(N - I)} 1 ' 2 n 2 z \N -2:0ft) (58) 

we are able to express the probability in terms of integrals over solid angles: 

a n J a n n z n z (n ■ n ) 

(59) 



P2 



l >] ' d 2 n / d 2 n' n z ri z (n ■ n') N ~ l 



12 



We suggest the following procedure to calculate these integrals. One first integrates over ft' for 
a fixed n , using spherical coordinates relative to the 'vertical' axis directed along n : the polar 
angle 0' is then the angle between n' and n so that one has simply n ■ n' = cos 9' . The 
integral over 6' and over the azimuthal angle (/)' can be performed, giving a result involving 
only n z . The remaining integral over n is performed with the spherical coordinates of vertical 
axis z . This leads to 

The ratio P2/P1 is therefore different from the naive (and wrong!) prediction (|55|) . For 
N = 2 one finds P2/P1 = 1 so that the second atom is surely in m = if the first atom was 
detected in m = . As the two atoms were initially in the state with total angular momentum 
zero, this result could be expected from the expression (|19|) of the two-particle spin state. In 
the limit of large ./V we find that once the first atom has been detected in the m = channel, 
the probability P2/P1 for detecting the second atom in the same channel m = is 3/5 . This 
somehow counter-intuitive result shows that the successive detection probabilities are strongly 
correlated in the case of the spin state (fl9|) . 

The exact calculation of the ratio 



Pk+i 



d 2 n I d 2 n n k+1 n' k+1 (n ■ n') N ' {k+1) 

(61) 



' d 2 n d 2 n' n k n' k (n-n') N ~ k 



gets more difficult when k increases. The large N limit for a fixed k is easier to obtain, para- 
phrasing the reasoning leading to Eq. (|13[) : in the integral over n' the function (n ■ n') Ar_ ( fc + 1 ) 
is extremely peaked around n' = ft so that we can replace the slowly varying function n' k+l 
by n k+1 . This leads to 

Hm »±i = a*±i. (62) 

n^+oq p k 2k + 3 

We now give the reasoning in the symmetry breaking point of view (|36|) . which assumes that 
a single experimental realization of the condensate corresponds to a Hartree-Fock state \N : 0,^) 
with the direction n being an unpredictable random variable with uniform distribution over 
the sphere. If the system is initially in the spin state \N : 0^) there is no correlation between 
the spins, and the probability of having k detections in the channel m = is simply (n 2 ) k . 
One has to average over the unknown direction n to obtain 

One recovers in an easy calculation the large N limit of the exact result, Eg. 1)62(1 ! We note 
that the result (|63|) is much larger than the naive (and wrong) result l|55J) as soon as k 3> 1 . 
In the limiting case where all the atoms of the gas have been detected, that is k = N , we 
expect to see a difference between the symmetry-breaking prediction and the exact result. The 
probability p^ of finding all the N spins with vanishing angular momentum along z is given 
for the exact ground state by 

^ = |(iV:O e -J^o)| 2 = F ^ T (64) 
so that the symmetry breaking prediction p s ^ is too small by a factor of two. 
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3 The ground state of a ID attractive Bose gas 

We consider in this section a model of a one-dimension Bose gas with attractive interactions in 
the absence of confining potential. This model is not so unrealistic at it may appear at a first 
glance. Experimentally one can start with a three-dimensional Bose-Einstein condensate with 
effective attractive interactions, that is with a negative scattering length, as it is the case for 
lithium [3j or rubidium 85 |2E]. The condensate is then subject to a strong harmonic confinement 
in the x — y plane, with a quantum of oscillation along x and y larger than the absolute value 
of the mean interaction energy per particle in the gas. In this way the motion of the condensate 
atoms is frozen transversally in the ground state of the harmonic trap in the x — y plane. Along 
the z direction one slowly reduces the trap strength so that the gas becomes almost free. 

Such a situation is interesting physically as it gives rise to the formation of 'bright' solitons 
well known in optics but not yet observed with atoms. Moreover we have found exact TV— body 
eigenstates for the model of a one-dimension Bose gas with a 5 attractive interaction potential 
in free space. We use these exact solutions to test the translational symmetry breaking Hartree- 
Fock approximation. 

3.1 A model for the interaction potential 

By analogy with the three-dimensional model interaction potential (|14|) we model the interaction 
in the one-dimension Bose gas by a contact potential: 



V(z 1 - z 2 ) = gid${zi - 



(65) 



The validity of such a modelization, and the value of the coupling constant gid in terms of the 
three-dimensional scattering length a are discussed in |27j . 

The 5 potential has the advantage in ID of leading to a mathematically well defined scat- 
tering problem, so that no regularization operator is required, contrary to the 3D case. In a 
practical solution of the N— body problem the use of a 5 potential amounts to imposing some 
boundary conditions on an otherwise interaction free Schrodinger equation |28[ I29j. Consider 
indeed an eigenstate of the N— body Hamiltonian with energy E : 



E^(zi,...,z N ) 



2m 



2 N 



£^+ 9ld £ 6(2 

l<i<j<N 



Zj) 



8=1 



. . .,Z N ). 



(66) 



The N— body wavefunction ^ is bosonic: it is symmetric with respect to any permutation of 
the N coordinates Z\, . . . , zjy so it is sufficient to determine it on the fundamental domain 



D = {(z\, . . . , zn) such that Z\ < . . . < zn}. 



(67) 



Inside the domain D the positions of the particles are different so that the 5 terms can be 
omitted in (|66|) : 

-2 N 



2m 



i=l 



* = 



(68) 



which is Schrodinger's equation for a non-interacting gas. The 5 terms are responsible for a 
discontinuity of the spatial derivatives of across the boundary of D . Using again the bosonic 
symmetry of ^ one is able to relate a left derivative (that is a derivative evaluated from the 
external side of D ) to a right derivative (that is a derivative evaluated from the internal side 
of D ) so that one obtains the conditions: 



3 Z , f 



=0+ 



h 2 



z 3 + l= z j 



(69) 
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The differential equation Q68|) over domain D , the conditions (|69f) on the boundary of D 
and the condition that |^| remains finite even for particles going to ±00 define the eigenstates 
of our model. 

3.2 Ground state of the one-dimension attractive Bose gas 

It turns out that in the previously exposed model one can calculate exactly eigenenergies and 
eigenstates of the Hamiltonian for N particles using the Bethe ansatz |28l l2"9 in the fundamental 
domain D : 

I N \ 

fy(zi,... ,zn) =N 22 ^(o-)exp I i22zjk a (j)\ for (z\, . . . , zn) G D (70) 
°es N \ j=i J 

where the kj 's are real arbitrary wavevectors and where the sum is taken over all permutations 
a of N objects. This form is a superposition of plane waves in the domain D and clearly 
solves l|5H|) with the energy 

h 2 N 

E = —Y / kl (71) 

The amplitudes A(a) of each permutation in the sum are adjusted to satisfy the boundary 
conditions (|6T)|) : 

m- n (> + f^t)- < 72 > 

i<j<i<N \ Mj) MO/ 

In the case > one can show that the family of eigenstates obtained from the Bethe ansatz 
is complete j^H] • The proof of completeness fails for the attractive case g\d < . 

We are considering here the less studied attractive case g±d < . Intuitively, the fact that 
this family of eigenstates is no longer complete is not surprising as the Bethe ansatz leads to 
positive eigenenergies only, whereas bound states with negative energies are expected to exist 
for attractive interactions. The way out is to extend the Bethe ansatz by adding an imaginary 
part Kj to the wavevectors: 

kj — > Kj = kj + iKj (73) 

where both kj and Kj are real. After this substitution in the Bethe form (|70[) the N— body 
wavefunction still satisfies Ea.(|68|) with an energy 

h 2 N 

The boundary conditions (|69j) are equally satisfied provided that one performs the substitution 
(|73|) in Ea.(|72j). What remains to be verified is that the N— body wavefunction does not diverge 
exponentially when one or several particles are going to infinity! For the term of the generalized 
Bethe ansatz associated to the permutation a , this requires that either 

A(a) = (75) 

or the real part of the argument of the exponential in the Bethe ansatz should not run off to 
+00 over the domain D |H0| : 

TV 

ZjK a (j) > —00 over D. (76) 
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The second condition (|76[) can be rewritten in a more explicit way, by taking z\ (varying from 
— oo to +00 ) and Z2 = Z2 — zi, . . . , Zn = zn — zn-i (varying from to +00 ) as independent 
variables. Equation l|76|l is then equivalent to the set of conditions 

N 

$>i = (77) 

N 

Z for l = 2,...,N. (78) 

3=1 

Once these conditions are satisfied one can hope to have a vanishing imaginary part of the energy 
E in Eq.dZU)! 

We have found a general form for the Kj satisfying for all a 's one of the two conditions 
(|75I76|) . We discuss here the set of Kj corresponding to the absolute ground state of the gas, 
the case of excited states is considered in fcj.51 First it is found that all the real parts kj are 
vanishing, which indeed is an efficient way to minimize the energy E as the kj 's contribute as 
squares with a positive coefficient! Second the imaginary parts Kj are given by 

^ = ^^[2j-{N + l)] for j = l,...,N. (79) 
The corresponding ground state energy was already known |31j : 

E (N) = -±^N(N 2 -l). (80) 

Let us check that the choice (|79*)l indeed leads to a N— body wavefunction not diverging at 
infinity. 

A first way to proceed is to show that conditions ([75 1) or (|76|) are satisfied. First one remarks 
that a permutation a different from identity will have a vanishing coefficient A(a) . The 
reciprocal a~ 1 of such a permutation cannot indeed conserve the numerical ascending order of 
the integers 1, . . . , N so that there exists an integer n satisfying 

a- 1 (n+ 1) < a- l (n). (81) 

If one considers in A(a) the factor of indices j < I with 

j = (j- 1 ^ + l) and I = cr _1 (n) (82) 

one gets a vanishing factor in Eg. 1)72(1 : 

1 + ; m9id/h2 . = (ss) 

according to Eq.(|79jl. leading to A(a) = . Second one can check by direct substitution that 
conditions (|77l78f) are satisfied, using the identity 

N 

J2[2j-(N + l)] = (N + l-l)(l-l). (84) 
3=1 

A second way to show that the proposed ground state wavefunction cannot blow up at 
infinity is to realize that, up to a normalization factor N , it can be rewritten as 



^o(zi, - ■ ■ ,z N ) =AAexp 



l<i<j<N 



(85) 
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for all values of the coordinates Zj 's, not restricting to the fundamental domain D . This form, 
already known in the literature, coincides with the generalized Bethe ansatz over the fundamental 
domain D and is obviously symmetric with respect to any permutation of particles. As the 
argument inside the exponential in (jHSJ) is negative, the expression (j85j) is clearly bounded from 
above. To determine the normalization factor J\f we enclose the gas in a fictitious box of size 
L tending to +00 ; this leads to [32] : 



{N-l)\ fm\g ld 



NL 



N-l 



(86) 



3.3 Hartree-Fock approximation for the ground state 

To what extent can we recover the results of the previous subsection using a Hartree-Fock ansatz 
\N : (f>) for the ground state wavefunction? We calculate the mean energy of the Hartree-Fock 
ansatz, and we obtain the following energy functional of <j> , sum of kinetic energy and mean 
field interaction energy: 



E 



hfr 



j dz 


'h 2 


dc/) 


2 


dz 


+ 




2m 



(87) 



This functional has to be minimized over cf) with the constraint that cf) is normalized to unity. 
A pure dimensional analysis gives already the main feature of the result. Let us rescale the 
coordinate z with some length £ : 



(j>{z) = -^(x = z /0 

such that the kinetic energy term and the mean field energy term have equal coefficients: 

h 2 (N-l)\gu\ 



m^ 2 



so that the energy functional becomes 



m 9u 



dx 



N(N- l) 2 e[-0,V*] 



(88) 

(89) 

(90) 
(91) 



We have already obtained the minimal Hartree-Fock energy up to a numerical factor. 

The minimization of the energy functional e gives access to this numerical factor. One 
can restrict to a real wavefunction ip(x) since a x— dependent phase for a fixed modulus 
immediately increases the kinetic energy without decreasing the interaction energy. The fact that 
the functional e is stationary around the minimal energy solution ip leads to the adimensioned 
Gross-Pitaevskii equation: 

~ijj"(x)-if) 3 (x) = u^(x) (92) 

where the adimensioned chemical potential v is a Lagrange multiplier ensuring that ip remains 
normalized to one in the variation. After multiplication of (|92j) by ip'(x) one can integrate over 
x to obtain: 

- \^' 2 {x) - \^{x) - \v^\x) 







(93) 
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where the integration constant has been taken equal to zero as tp vanishes at x = ±00 . Note 
that the quantity v is then clearly negative. The problem is reduced to a quadrature and one 
finds: 

^x) = 2/25. (94) 

cosh[-^/2|f|(a; — xo)] 

where xq is arbitrary. As the integral over all the real axis y of 1/ cosh 2 (y) is equal to two, 
we obtain 

u = -\. (95) 

We have recovered the well-known solitonic solution for the one-dimension non-linear Schrodin- 
ger equation: 

2£V 2 cosh[(z - z )/{2Q\ 

parameterized by the arbitrary position zq of the soliton. Redimensioning the quantity u gives 
the chemical potential of the soliton [22] : 

^ = -\^{N-lf. (97) 
The final result for the Hartree-Fock minimal energy is of course independent of zq : 



Ef{N) = -^^N{N-lf. (98) 



The deviation of the Hartree-Fock minimal energy (|98|1 from the exact result ()8U(I is a fraction 
1/N of the energy and is small in the large N limit. Can we understand why the validity 
condition of the Hartree-Fock ansatz is simply N S> 1 ? The result is not intuitive as large 
values of N lead to small spatial widths £ of the soliton, and therefore to high linear densities 
of particles, where one may expect to have a strongly interacting regime not well described by 
mean field. The paradox can be removed in the following way: the Hartree-Fock ansatz for 
the N— body wavefunction does not contain any correlation between the particles, as it is a 
factorized state vector, so it is an acceptable approximation only if the interaction potential is 
weak enough, that is if the interaction potential can be treated in the Born approximation for 
the relevant relative momenta of the particles. An exact calculation of the scattering amplitude 
of two particles with initial relative wavevector k and interacting with the model potential (165(1 
shows that the Born approximation is applicable provided that k is high enough: 



h 2 k 



mg ld 



> 1 (99) 



in contrast to the three-dimensional case. Taking the estimate k — l/£ for the typical relative 
momentum of particles we find that condition (|99|) reduces to N 3> 1 ■ Actually the fact that 
the weakly interacting regime in one-dimension corresponds to a high density regime of the gas 
is a well established fact also for effective repulsive interactions g\<i > |28j . 

There is however a notable difference of translational properties of the exact TV— body 
ground state (|85|) and of the Hartree-Fock ansatz. Whereas the exact ground state is invariant 
by a global translation of the positions of the particles, as it should be, the Hartree-Fock ansatz 
leads to condensate wavefunctions 4> ZQ localized within the length £ around some arbitrary 
point zq , see Eq.(|96|). The Hartree-Fock ansatz \N : <j) zo ) therefore breaks the translational 
symmetry of the system. 

Breaking a symmetry of the system costs energy, and this can be checked for the present 
translational symmetry breaking. As the center of mass coordinate and momentum Z, P of 
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the N particles are decoupled from the relative coordinates of the particles we can write the 
total energy of the gas as the sum of the kinetic energy of the center of mass and of an 'internal' 
energy including the kinetic energy of the relative motion of the particles and the interaction 
energy. Whereas the exact ground state wavefunction has a vanishing center of mass kinetic 



energy, the symmetry breaking state \N 



J z / 



E 



hf 



(N:<f>, 



contains a center of mass kinetic energy: 



(100) 



where mN is the total mass of the gas and P = p\ + . . .+pn is the total momentum operator 
of the gas. Expanding the square of P , and using the fact that the soliton wavefunction <j> has 
a vanishing mean momentum we obtain 

„2 



E 



hi 



(0 



V 

Zol 2m' 
1 mg\ d 



(N - If 



(101) 
(102) 



24 h 2 

We see that E~r o m accounts for half the energy difference between the exact ground state 
energy (p")|) and the Hartree-Fock energy (|9*8)) . 
The restored symmetry ansatz of Eq. Q : 



1% 



vhfs 



r N 



L/2 
-L/2 



dz \N : (j) z 



(103) 



has the advantage of having a vanishing total momentum and will have a lower energy than 
the symmetry breaking ansatz. We can try to calculate its energy. First we have to evaluate 
overlap integrals of 4> ZQ and 4> z ^ , with zq ^ z' : these integrals can be brought in the form of 
integrals of rational functions using a change of variable z = 2£ exp u ; these integrals over u 
can be calculated exactly. We are left with an integral over zq = x/(2£) : 

, N 

dx A(x) 



E^(N) = E?(N) 



hf/ 



A(x) 



x 



sinh x 



dx 



N 



3 + cosh 2x 



sinh a; 

2 sinh 2x 



+ 



6 



x sinh x 



cosh x 



sinhx 



(104) 



(105) 



2 sinh x 

This ratio of integrals can be evaluated in the large N limit, using Laplace's method. The 
function x/sinh(x) raised to the large power N is very peaked around x = with a width 
scaling as 1/ VJSf so that one can expand in powers of x = y / yN : 



N 



sinh x , 
A(x) 



exp 



i+§£+o<w 



We are left with Gaussian integrals that can be performed exactly, 
results to the exact ground state energy: 

2 



E™(N) 



E^(N) 



E (N) 
Eq{N) 



1 



N 
4 

5N 



+ 0(N- 2 ) 
+ 0(N- 2 



(106) 
(107) 

We compare the various 

(108) 
(109) 



We see that the symmetrized Hartree-Fock prescription @ has a lower energy than the symmetry 
breaking prescription (JJJ) but the relative error is of the same order 0(1/N) . 
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3.4 Physical advantage of the symmetry breaking description 

We now raise the question: is there a Bose-Einstein condensate in the one-dimension unconfined 
Bose gas with attractive interaction? To make things simple we assume that the gas is at zero 
temperature so that the N— particle wavefunction is known exactly, see Eq.(|85|). 

We start with a reasoning in terms of the one-body density operator (even if we know from 
the previous physical example that this may be dangerous, see ij2.4|) . Paraphrasing the usual 
definition of a Bose-Einstein condensate in three dimensional free space we put the one-dimension 
gas in a fictitious box of size L and we calculate the mean number of particles no in the plane 
wave with vanishing momentum p = in the limit L —* +oo . 

The calculation with the exact ground state wavefunction can be done, it is too involved to 
be reproduced here. One finds that no is going to zero as 1/L : 

C(N) 2h 2 

n o - —r-^—, r- ( 110 ) 

L m\gu\ 



The factor C(N) is given by 



k(N + l-k)-^(N + 1) 



1111 



and converges to 7r 2 /2 in the large N limit. There can be therefore no macroscopic population 
in the p = momentum state in the large L limit. One may then be tempted to conclude 
that there is no Bose-Einstein condensate, even at zero temperature, in the free one-dimension 
Bose gas with attractive contact interactions. However we have learned that a reasoning based 
on the one-body density matrix may miss crucial correlations between the particles, and that 
the symmetry breaking point of view may be illuminating in this respect. 

The translational symmetry breaking point of view approximates the state of the gas by the 
N— body density operator: 

P sh = lim / -^|iV:^ o )(iV:0 Zo |. (112) 



In the large N limit we expect this prescription to be valid for few-body observables. Of course 
for a N— body observable such as the kinetic energy of the center of mass of the gas, the results 
will be different: the kinetic energy Ea. HlU2|) for the symmetry-breaking point of view differs 
from the exact vanishing value. 

Let us test this expectation by calculating in the Hartree-Fock approximation the mean 
number of particles in the plane wave {z\k) = exp(ikz)/ X 1 / 2 . Using the following action of the 
annihilation operator of a particle with wavevector k on the Hartree-Fock state: 

a k \N : <p zo ) = N x l 2 {k\<j> Z0 )\N - 1 : ^ ) (113) 

we obtain 

nf = NM^)? . (114) 

The momentum distribution of the particles in the gas in this approximation is simply propor- 
tional to the momentum distribution of a single particle in the solitonic wavefunction <j) Zo ! It 
turns out that the Fourier transform of the 1/ cosh function can be calculated exactly, and it 
is also a 1/ cosh function. We finally obtain: 

n£ f ~-— r ^ 115) 

Lm\g ld cosh 2 (7r£;£ 
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where £ is the typical soliton size given in Ea.(|89[). For k = one recovers the large N limit 
of the exact result (|11U|) . 

In more physical terms, one can imagine from Eg, 1)112(1 that a given experimental realization 
of the Bose gas corresponds to a condensate of ./V particles in the solitonic wavefunction (|96|) . 
with a central position zq being a random variable varying in an unpredictable way for any 
new realization of the experiment. There is therefore a Bose-Einstein condensate in the one- 
dimension attractive Bose gas! 

An illustrative gedanken experiment would be to measure the positions along z of all the 
particles of the gas. In the symmetry breaking point of view the positions zx,...,zn obtained in 
a single measurement are randomly distributed according to the density |(^ | (z) = \4>o( z — z o)\ 2 
where zq varies from shot to shot as the relative phase of the two condensates did in the 
MIT interference experiment. As we know the exact ground state (|85|) we also know the exact 
A— body distribution function, l^o^i, • • • > z n)\ 2 ■ This is however not so easy to use! 

So we suggest instead to consider the mean spatial density of the particles knowing that the 
center of mass of the cloud has a position Z . In the exact formalism this gives after lengthy 
calculations: 

p{z\Z) = J W.. J ' dz N \M*i,---,ZN)\ 2 (^5(z- Zj ))L5(z-^J2 z n) ( 116 ) 



N ^ (A -2)! Nl . lVfc . ; 



i) v 



A — 1 £ 



p sh (z 



where £ is the A— dependent length of the soliton (|89|). the integrals are taken in the range 
[—L/2,L/2] and L —* +oo ; the factor L , compensating the one in the normalization factor 
of ^ , ensures that the integral of p{z\Z) over z is equal to N . 

In the symmetry breaking point of view the definition of p(z\Z) is similar to Ea. H116|) ; the 
factor L cancels with the 1/L factor of Eq. (|112j) . This leads to 

Z) = JdzoJd Zl ..Jdz N (jl\c/> Zo (z k )\ 2 ^\^S(z-z^ 

= N ] d Zl ... j dz N (^\Mzk)\^s{z-z + Zl -^Y.^ 

where we have made the change of variables Zk — > z^ + zq (which allows to integrate over 
zq ), we recall that c^o is the solitonic wavefunction centered in zq = and we have replaced 
the sum over the indistinguishable particles j by N times the contribution of particle j = 
1 . The multiple integral over the positions Z\,...,Zn can be turned into a single integral 
over a wavevector q by using the identity 5(X) = J dq/(2ir) exp(iqX) , allowing a numerical 
calculation of p sh (z\Z) . 

Does the approximate result (|117j) get close to the exact result for large A ? We compare 
numerically in figurenthe exact density p(z\Z) to the symmetry breaking mean-field prediction 
p sb (z\Z) : modestly large values of A give already good agreement between the two densities. 
This validates the symmetry breaking approach for the considered gedanken experiment. 

What happens in the large A limit? In Eq. (|117|) each variable z^ explores an interval of 
size ~ £ so that the quantity (z\ + . . . + zn)/N has a standard deviation ~ £/V~N much 
smaller than £ and can be neglected as compared to z\ inside the 5 distribution. This leads 
to _ 

p sb (z\Z) ~ A|<^ 0=z (*)| 2 for VA»1. (118) 

Numerical calculation of p sh (z\Z) shows that Ea. (|118|) is a good approximation over the range 
\z — Z\ ~ £ for A = 10 already! 
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Figure 1: For the ground state of the one-dimension attractive Bose gas, position dependence of the 
mean density of particles knowing that the center-of-mass of the gas is in Z = . Solid line: exact 
result p(z\Z = 0) . Dashed line: mean- field approximation p s] °(z\Z = 0) . The position z is expressed 
in units of the 'soliton' radius £ where £ is given in Ea. (j89fl . and the linear density in units of iV/£ . 
The number of particles is (a) N = 10 and (b) N = 45 . 



3.5 Excited states of the one-dimension Bose gas 

We briefly show how to calculate excited states of the gas which can contain clusters or 'lumps' 
of bound particles, and that are therefore not contained in the usual Bethe ansatz (|7U|). These 
excited states are obtained with the generalized Bethe ansatz involving complex wavevectors 
Kj as defined in Ea.(|73j). 

The first step is to group the atoms in an arbitrary number R of lumps ( 1 < R < N ), 
the lump number r containing n r particles. More precisely, the positions of the particles are 
in the fundamental domain D , that is they are in ascending order z\ < . . . < zn ; we put the 
particles 1, . . . , n% in the first lump, the particles n\ + 1, . . . , n\ + ni in the second lump, etc, 
the particles N — ur + 1, . . . , N in the last, R th lump. It is then more convenient to reindex 
the particles: each atom is now identified by its lump number r (from 1 to R ) and its rank 
inside the lump i r (from 1 to n r ) 

The complex wavevectors are then taken as 



K 



r,i r 



m\gi d \ 
2h 2 



[2i r - (n r + 1)]. 



(119) 



as examplified in figure |21 This choice has a simple physical interpretation. The imaginary parts 
of the Kj 's inside lump r with n r particles are exactly the Kj 's defining the ground state of 
a gas of n r atoms. The real part of the Kj 's inside lump r is the same for the n r particles 
of the lump and corresponds to a global motion of the lump with a total lump momentum of 
n r hk r : this is nothing but an excitation of the center of mass motion of the particles in the 
considered lump. The corresponding N— body wavefunction is then simply a coherent ensemble 
of R quantum solitons with different momenta Hk r . The eigenenergy associated to the choice 
(|119|) can be calculated from formula (|74|): 



2 ,X 



E 



E (n r ) + n, 



2m 



(120) 



where the ground state energy Eq is given in (|8()j) as function of the number of particles. We 
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Figure 2: We plot the momenta Kj 
particles with 7 lumps. 



kj + iKj for an arbitrary excited state of N = 21 



therefore find that the R quantum solitons do not interact: this is the quantum analog of the 
known fact that classical field solitons can cross without interaction. 

What remains to be checked is that the proposed N— body wavefunction does not explode 
at infinity. The proof proceeds along the same lines as for the ground state. One first shows 
that only the permutations a such that <t _1 conserves the ascending order of particle labels 
inside each lump can have a non-vanishing factor A(a) . That is a has to satisfy 



a 



~Ht,1) < ... < a l (r,i r ] 



< 



< a~ 



[r, ?v 



for all 



l,...,R. 



:i2D 



Then the conditions H77I78|) are checked to be satisfied by direct substitution. 

It is interesting to compare the exact spectrum (j!20j) to mean field predictions. We restrict 
here [3D to a simple calculation using the Bogoliubov approach to calculate the excitation spec- 
trum of the classical soliton. This amounts to linearizing the time dependent Gross-Pitaevskii 
equation 



ihdt4>(z, t) 



2m 



d 2 z <f>(z, t) + g ld {N - l)\<J)(z, t)\ 2 (j)(z, t) - n<t>(z, t) 



(122) 



around the stationary solution (|96|) and to looking for corresponding eigenmodes. One finds that 
the eigenmode (u(z),v(z)) with energy e solves the usual Bogoliubov-de Gennes equations: 

^ d 2 



€u(z) 



-ev(z) 



A* + 20n(]V-l)|* 
Zm dz A 



2m dz 2 



+ 2g ld (N-l)\<j> z 



I'- 



ll 



u{z)+g{N -l)4> 2 zo v(z) (123) 
v{z)+g{N-l)4>fu(z). (124) 



The Bogoliubov modes have long been known in the context of optical solitons |35j . Apart 
from two zero energy modes corresponding to the £7(1) symmetry and translational symmetry 
Goldstone modes, the eigenmodes behave as traveling waves exp(ikz) far away from the center 
of the soliton. The eigenenergy is straightforward to obtain: far from the soliton we can set 
<t>za — in (|123fl and we get the dispersion relation 



h 2 k 2 

2m 



(125) 
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corresponding to a gap \fx\ in the excitation spectrum. In the large N limit the Hartree-Fock 
chemical potential /j, of (|97|) is very close to the exact chemical potential Eq(N) — Eq(N — 1) . 
In this way the elementary excitation energy is very close to the energy difference between 
the absolute ground state of the gas and the excited state with two lumps, one lump with N — 1 
particles at rest and the other lump with one particle with momentum fik ! The Bogoliubov 
approach is however unable to predict bound states between excitations, corresponding to excited 
lumps with two or more than two particles. 

4 Conclusion 

We have studied in great details two exactly solvable models of the interacting Bose gas. The 
first model involves N interacting spin one particles assumed to be in the ground state of a 
cubic box, the second model is a one-dimension Bose gas of spinless particles interacting with 
an attractive contact potential in the absence of external confinement. 

Both models share the property of not leading to the formation of a Bose-Einstein condensate, 
even at zero temperature, in the strict Penrose and Onsager sense: the one-body density operator 
has three macroscopic eigenvalues for the spins, and has no macroscopic eigenvalue for the 
attractive Bose gas. 

We have shown that the usual Hartree-Fock approach, trying to approximate the ground 
state of the gas by a Hartree-Fock state with all the N particles in the same single particle 
state (ft , can bring a considerable physical understanding of these model systems. A continuous 
family of (ft 's is found to minimize the mean energy, each member of the family having the 
same energy and corresponding to a state with broken rotational symmetry (for the spins) or 
translational symmetry (for the one-dimension gas). This allows to imagine that any particular 
experimental realization of the gas is a condensate in the broken symmetry state (ft , with (ft 
varying in an unpredictable way for any new realization of the experiment. E.g. one imagines 
that any particular realization of the one-dimension attractive Bose gas is a condensate in a 
soliton wavefunction. 

We have successfully tested the Hartree-Fock approximation against the exact results, in the 
thermodynamical limit for the spins and in the large N limit for the one-dimension Bose gas. 
In particular the prediction of the issues of gedanken experiments performed on the gas was 
found to be much more transparent in the symmetry-breaking point of view than with the exact 
N— body ground state wavefunction. 

An interesting aspect has not been discussed in the present article. It involves the compared 
robustness of the exact ground state and the Hartree-Fock states against decoherence induced e.g. 
by particle losses due to three body collisions jSHI- While the Hartree-Fock states \N : (ft) remain 
Hartree-Fock states under removal of one or several particles and are therefore robust, the exact 
ground state, after having experienced several particle losses, becomes a density operator that 
may actually look like a statistical mixture of Hartree-Fock states, these losses being formally 
equivalent to measurements performed on the system |10j . Such an aspect is important in the 
context of an experimental attempt to distinguish between the exact and the symmetry breaking 
descriptions. 

We acknowledge useful discussions with Gordon Baym, Jean Dalibard, Anthony Leggett, 
Maxim Olshanii and Sorin Paraoanu. We are grateful to Alice Sinatra for a critical reading of 
the manuscript. 
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